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Abstract 



We define the Schrodinger equation with focusing, cubic nonlinearity on one-vertex graphs. We prove 
Oh global well-posedness in the energy domain and conservation laws for some self-adjoint boundary condi- 

tions at the vertex, i.e. KirchhofF boundary condition and the so called S and S' boundary conditions. 
Moreover, in the same setting we study the collision of a fast solitary wave with the vertex and we show 
that it splits in reflected and transmitted components. The outgoing waves preserve a soliton charac- 
G ter over a time which depends on the logarithm of the velocity of the ingoing solitary wave. Over the 

same timescale the reflection and transmission coefficients of the outgoing waves coincide with the cor- 
f>r^ responding coefficients of the linear problem. In the analysis of the problem we follow ideas borrowed 

^ from the seminal paper [17j about scattering of fast solitons by a delta interaction on the line, by Holmer, 

Marzuola and Zworski; the present paper represents an extension of their work to the case of graphs and, 
^0 as a byproduct, it shows how to extend the analysis of soliton scattering by other point interactions on 

the line, interpreted as a degenerate graph. 

\l Keywords: quantum graphs, non-linear Schrodinger equation, solitary waves. 
§ MSC 2010: 35Q55, 81Q35, 37K40. 

1 Introduction 

X 

In the present paper we study the nonlinear wave propagation on graphs. As far as we know the subject 
of nonlinear Schrodinger evolution on graphs is at its beginnings. An extensive literature on the behaviour 
of linear wave and Schrodinger equations on graphs exists ([22l [23l [2T1 [TJ [8], and references therein) and a 
certain activity concerning the so called discrete nonlinear Schrodinger equation (DNLSE) in chains with 
edges and graphs inserted ("decorations", interpreted as defects in the chain), both from the physical and 
the numerical point of view (see for example [201 E]). To our knowledge, however, there are only very few 
papers in which nonlinear Schrodinger evolution on graphs has been introduced and studied (see |6l |26] ) . 
In the first paper [B] (and similar ones quoted therein) NLS on graphs emerges in some models of quantum 
field theory on "bulks" ; the second recent paper [26j addresses from a physical point of view some general 
questions related to the ones here studied and is briefly commented in the conclusions of the present paper. 
Besides several results on existence and stability of stationary states for the NLS on star graphs can be 
found in [2j. 
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In any case the study of nonlinear propagation in ramified structures could be of relevance in several 
branches of pure and applied science, from condensed matter physics to nonlinear fiber optics, hydrodynamics 
and fluid transport (a non-traditional example is blood flow in veins and arterias), and finally neural networks 
(see for example the study of reaction- diffusion type FitzHugh-Nagumo-Rall equations on networks in jl2] . 
and references therein). 

In all these examples there is a strong dependence on the modelization. The nonlinear Schrodinger equa- 
tion with cubic nonlinearity is especially suitable to describe nonlinear electromagnetic pulse propagation in 
optical fibers and, under the name of Gross-Pitaevskii equation, the dynamics of Bose-Einstein condensates. 
Better suited for other applications, for example the hydrodynamic flow, is the KdV equation or its relatives, 
not treated here. 

Of course real networks are not strictly one dimensional, and an abstract graph, which is just a set of 
copies of M"*" ("edges", or "branches") with functions living on it satisfying certain boundary conditions at 
lacks some of the geometric meaningful characteristics of a real network, such as thickness and curvature 
of the branch or orientation between edges. On the other hand, problems related to wave propagation on 
networks are far from being well understood also for the linear propagation (see for example [U [10] ) , and 
so we content ourselves with posing and analysing the nonlinear problem in the idealized and simplified 
framework of an abstract graph. 

We study here the special case of a star graph with three edges. A generalization to a star graph with 
n edges would be straightforward, but here our interest is in clarifying the main features of the evolution 
and the techniques involved in its analysis. A preliminary point and not a trivial issue is the definition of 
the dynamics. Let us recall that for a star graph Q, the linear Schrodinger dynamics is defined by giving a 
self-adjoint operator H on the product of n copies of L2(M+) (briefly L'^iG)), with a domain D{H) in which 
appears a linear condition involving the values at of the functions of the domain and of their derivatives. 
The admissible boundary conditions characterize the interaction at the vertex of the star graph. For the 
nonlinear problem, we establish the well-posedness of the dynamics in the case of a star graph with some 
distinguished boundary condition at the vertex, namely the free (or Kirchhoff), and the 6 and 6' boundary 
conditions (see section 2 for the relevant definitions). 

To clarify the exact problem we are faced to, the differential equation to be studied is of the form 

ij^^t = H^t - I'i'tl^'^t , t^O, (1.1) 
where the function ^ is, for a three-edge graph, a column vector 




that lies in the domain D{H) of the linear Hamiltonian H on the graph, so emboding the relevant boundary 
conditions. This is the abstract strong form of the equation, which is equivalent to a particular nonlinear 
coupled system of scalar equations. The coupling is not due to the nonlinearity, because of the definition 

/ l^iPV'i \ 

but to the boundary condition at the vertex. For example, in the simple case of a Kirchhoff boundary 
condition, the coupling between the edges is given by 

^ G L\g) s.t. t/ji G H\R+), ^i(O) = ^2(0) = ^3(0), V'i(O) + V'aCO) + V'sCO) = . 

The 5 or 5' boundary conditions allow a coupling between the values of the function \I' and the values of 
their derivatives at the origin, but in principle the nature of the problem is unaltered. 
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In the present paper, for several reasons, we prefer to write the dynamics in weak form, which is the following 



$i = e-'"^^o + i [ e-'^'^^-'^l^sl^'^sds, t^O; (1.2) 
Jo 

here the ^ belongs to the form domain D^Sun), where Sun is the quadratic form of the linear Hamiltonian 
H. We interpret, according to the use, the form domain as the finite energy space. An adaptation of 



the methods in |3] gives the local well-posedness of the equation (1.2) for every initial data in 
Moreover, charge and energy conservation laws hold true for such weak solutions, and as a consequence, the 
NLS on graph admits global solutions. The generalization of the well-posedness result to the case of general 
self-adjoint boundary conditions will be treated elsewhere. 

Apart from well-posedness, the main goal of this paper is to provide information on the interaction 
between a solitary wave and the boundary condition at the vertex. As it is well-known, the NLS on the line 
admits a family of solitary non dispersive solutions, rigidly translating with a fixed velocity. A rich family 
of solitary solutions ( "solitons" ) is given by the action of the Galilei group on the function 

(l){x) = \/2cosh"^x, xGM, 

or explicitly, 

2 

(pxo,v{x, t) = e*^^e"**'^re**(/.(a; - xq - vt) x£R,t£R,v£R. 

We show that, after the collision of a solitary wave with the vertex there exists a timescale during which 
the dynamics can be described as the scattering of three split solitary waves, one reflected on the same branch 
where the originary soliton was running asymptotically in the past, and two transmitted solitary waves on 
the other branches. On the same timescale, the amplitudes of the reflected and transmitted solitary waves 
are given by the scattering matrix of the linear dynamics on the graph. The soliton-like character persists 
over time intervals that depend on the velocity of the impinging original soliton: the faster is the original 
soliton, the (logarithmically in the velocity v) longer is the survival time of the solitary wave behaviour on 
every branch of the graph. The non-trivial point is that the persistence time of solitary behaviour after 
collision with the vertex is much longer, for fast solitons, than the time over which it is reasonable to 
approximate the nonlinear dynamics with the linear one. The same timescale of the order Inv of persistence 
of solitary behaviour appears in the paper [T] where the collision of two solitary waves with an underlying 
smooth potential is studied, and in the paper by Holmer, Marzuola and Zworski [17J on the fast NLS-soliton 
scattering by a delta potential on the line, which is the main source of inspiration for our result and for the 
techniques employed in the present paper. 

We give now an outline of our result and its proof. 
The initial data are of the following form 

(x(a:)e~*5^(/>(a; — xq) \ 
j , xo ^ v^-\ < (5 < 1, (1.3) 

where x is a cut off function, that is x ^ C°°(M^), x = 1 in (2,+oo) and x = in (0, 1). Apart from a 
small tail term truncated by the cutoff function, the first component is the initial condition of a free (i.e. 
without external potentials) NLS which on the line yields a solitary wave running with velocity v; the center 
Xq of the initial soliton is chosen far from the vertex. We are interested in the evolution of this initial 
condition. 

The dynamics can be divided into three phases. The pre-interaction phase, where the evolved initial 
condition is far from the vertex, and the undisturbed NLS evolution dominates. At the end of this phase, the 
solution enters the vertex zone, and differs (in norm) from the evolved solitary wave by an exponentially 
small error in the velocity v. The second phase is the interaction phase, in which a substantial fraction 
of the mass of the initial soliton has reached the vertex, and the linear dynamics dominates due to the 
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shortness of the interaction time, leaving the system at the end of this phase with three scattered waves, 
the amphtudes of which are given by the action of the scattering matrix of the associated hnear graph on 
the incoming sohtary wave. The size of the corresponding error is (again in L? norm) a suitable negative 
inverse power (t>~2'') of the velocity. In the phase one, the main technical tool consists in the accurate use, 
as fixed by |17j . of the Strichartz's estimates to control the deviations between the unperturbed NLS flow 
and the NLS flow on the graph. In the phase two, we need to compare the nonlinear evolution with the 
linear flow on the graph in the relevant time interval. Finally there is the post-interaction phase, where the 
free NLS dynamics dominates again; however, now the initial data are not exact solitary waves, but waves 
with soliton-like profiles and "wrong" amplitudes (due to the scattering process in the interaction phase). 

The true evolution is compared with a reference dynamics given by the superposition of the nonlinear 
evolution of the outgoing scattered profiles, and it turns out that the error is, in norm, of the order of 
an inverse power of velocity (depending on the size of the time interval of approximation). For a precise 
formulation one has to tackle the problem of representing the reference soliton dynamics to be compared 
with the true dynamics. This problem arises because one would like to use crucial and known properties 
of NLS on the line (such as existence of an infinite number of constants of motion), and various associated 
estimates, while on a star graph one has a NLS on halflines, jointly with boundary conditions. The problem 
occurs, of course, in each of the three phases in which the dynamics is decomposed. 

Our choice of reference dynamics is the following. We associate to every edge of the star graph a 
companion edge chosen between the other two, in such a way to have three fictitious lines; then we glue 
the soliton on every single edge with the right tail on the companion edge, respecting the free nonlinear 
dynamics. One of the main technical points in the analysis of the true dynamics is to have a control in the 
errors brought by this schematization. More precisely, let us define 



4't(xi,X2,X3) 



l'?(xi,X2,X3) 



re 



4 ig l^X2^tt^ 






{xi + xo - vt) ^ 

j(x2 - Xo + vt) 

I 



\ 

vt) 

'e~'^^^^e^^(p{xs - Xo + vt)/ 



4 *e*2^2gJt^(^2.2 _|_ xq 



;i.4) 



4 e ^i^^e^^cpixi - xq + vt) 







\ie *e*2^^e**i; 



(X3 + Xo - vt) 

Each of these vectors represents a soliton on the fictitious line given by an edge and its companion, multiplied 
by the scattering coefficients of the linear dynamics considered, Kirchhoff, 5 or 6' (and here left unspecified). 
Up to a small error, these functions represent outgoing waves at the end {t = t2) of the interaction phase, 
which is essentially a scattering process. Taking these as initial data for the free nonlinear dynamics on the 
pertinent fictitious line, we define their time evolution $j as given by 

ft 

i = 1,2,3, t^t2, 



Jt2 



where the Hj are the linear Hamiltonians that decouple the j + 2-branch from the others. With these 
premises, the main result of the paper is the following. 



Theorem 1.1. Let the unique, global solution to the Cauchy problem (1.1) with initial data (1.3). There 
exist r=K > and T* > such that for < T < T^: one has 



H - ^t\\Ll{g) 

for every time t in the interval t2 <t <t2 + Thiv , where C is a positive constant independent of t and v. 



^ Cv ^ 
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To be precise, the Hamiltonians in (1.1) to which the theorem refers have to be rescaled in order to give 



a nontrivial scattering matrix in the regime of high velocity (see section 4, in particular theorem 4.3). 

To get the previous result, Strichartz estimates do not suffice, and more direct properties coming from 
the integrable character of the cubic NLS are needed. In particular, thanks to the existence of an infinite 
number of integrals of motion, in [T7| a spatial localization property of the solution of NLS with smooth data 
is proven, with a polynomial (and not exponential) bound in time. This gives a control on the tails of the 
difference between the solution and the reference modified solitary dynamics. An analogous method applies 
in our case. Let us note that as a consequence of the previous result, we can give the scattered amplitudes 
in terms of the incoming amplitude and scattering coefficients of the linear dynamics, in the time range of 
applicability of the main theorem (see remark 4.4). 

We give a brief summary of the content of the various sections of the papers. In section 2 we give some 
generalities on linear dynamics on graphs, including Hamiltonians, their quadratic forms, resolvents and 
propagators. Moreover the essential Strichartz estimates are recalled. 

In section 3 local and global well posedness of nonlinear Schrodinger equations on star graphs is proved. 



In section 4 the main result (theorem 4.3) is introduced and stated. 

Section 5 is devoted to the proof of the result. In section 6 some final remarks are given and further possible 
developments are discussed. 

1.1 Setting and notations 

We consider a graph Q given by three infinite half lines attached to a common vertex. In order to study a 
quantum mechanical problem on Q, the natural Hilbert space is then L?'{Q) = L^(M"'') © L^(M"'") © Lp'iW^). 
We denote the elements of L'^{G) by capital greek letters, while functions in L^(]R"'") are denoted by lowercase 
greek letters. It is convenient to represent functions in L'^{Q) as column vectors of functions in L^(R~*"), 
namely 

/ i^i 

V V's 

The norm of L^-functions on Q is naturally defined by 



1*1 




Analogously, given 1 ^ r ^ oo, we define the space L^{G) as the set of functions on the graph whose 
components are elements of the space L'"(M+), and the norm is correspondingly defined by 
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ll*llL'-(e)= ' ^^r<oo, 11^-11^^ = sup IIV-jIIlo 

When a functional norm refers to a function defined on the graph, we omit the symbol G. Furthermore, 
from now on, when such a norm is L^, we drop the subscript, and simply write || • ||. Accordingly, we denote 
by (•, •) the scalar product in L^. 

As it is standard when dealing with Strichartz's estimates, we make use of spaces of functions that are 
measurable as functions of both time (on the interval [Ti,T2]) and space (on the graph). We denote such 
spaces by Lj*^^ rp^^U{Q), with indices 1 ^ r ^ cx), 1 < oo; we endow them with the norm 

/ rT2 \ 1/p 

ll^llz-f Jr(C\=\ I ,1^P<00, ll^llroo jr(c\= sup ll^'slLr- 

" "^[Ti,T2]^ ^'^) VVti / "^[Tl,T2]^ ie[Ti,T2] 

The extension of the definitions given above to the case p = c« r = oo is straightforward. 
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Besides, we need to introduce the spaces 
equipped with the norms 

3 3 

11*11^1(5) = 11*11/^2(5) = X] IIV'i|lH2(]K+). (1.5) 

i=l i=l 

The product of functions is defined componentwise, 

IV'ipVi 

2,T, - I l„/._|2„ 




'^^= I ^202 |, SO that l^'l"'^' EE I 1-021 V'2 

We denote by I the 3x3 identity matrix, while J is the 3x3 matrix whose elements are all equal to one. 
When an element of L?'{Q) evolves in time, we use in notation the subscript t: for instance, Sometimes 
we shall write ^{t) in order to highlight the dependence on time, or whenever such a notation is more 
understandable. 



2 Summary on linear dynamics on graphs 
2.1 Hamiltonians and quadratic forms 

Standard references about the linear Schrodinger equation on graphs are [71 13 [22l [231 [21] , to which we refer 
for more extensive treatments. Here we only give the definitions needed to have a self-contained exposition. 

We consider three Hamiltonian operators, denoted by Hp^ Hf , H^, (with a,/3 G M), and called, respec- 
tively, the Kirchhoff, the Dirac's delta, and the delta-prime Hamiltonian. These operators act as 




* ^ -i^'i (2.1) 



on some subspace of H'^{Q), to be defined by suitable boundary conditions at the vertex. 

Here and in the following subsection we collect some basic facts (see [21], [22], [23] [7]) on Hp, Hf, and 



The Kirchhoff Hamiltonian Hp acts on the domain 

V{Hf) := G H^iG) s.t. Vi(0) = MO) = M^), 0'i(O) + ^^2(0) + ^3(0) = 0}. (2.2) 



It is well-known, see |21j . that (2.2) and (2.1) define a self-adjoint Hamiltonian on L (Q). Boundary con- 
ditions in (2.2) are usually called Kirchhoff boundary conditions. We use the index F to remind that Hp 
reduces to the free Hamiltonian on the line for a degenerate graph composed of two half lines. 
The quadratic form £p associated to Hp is defined on the subspace 

v{£f) = {^€ H\g) s.t. Vi(o) = MO) = MO)} 

and reads 

3 r+00 
M*] = E / \i^'^i^)\'dx. 

i=l -^0 

The Dirac's delta Hamiltonian is defined on the domain 

V{H^) := G H\g) s.t. MO) = Mo) = Mo), Mo) + Mo) + ^^3(0) = aMO)} (2.3) 
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Again, is a self-adjoint opeator on L?'{Q) ([21]). It appears that Hf generalizes the ordinary Dirac's 
delta interaction with strength parameter a on the line, see, e.g. [5]. 
The quadratic form associated to Hf is defined on 



v{Ef) = {M/ G H\g) s.t. Vi(o) = V2(o) = V3(o)} 



and is given by 



" r+oo 

£m = Y. m^)\'dx+a\Mo)\'. 

The delta-prime Hamiltonian is defined on the domain 

P(i/f,) := G H\g) s.t. ^'i(O) = ^^(0) = ^^(0), V'i(O) + ^^2(0) + ^3(0) = M(0)} (2.4) 
Again, H^, is a self-adjoint opeator on L'^{Q) ([2T]). 

The quadratic form associated to H^, is defined on 'D{£^i) = H^{g) and is given by 



3 /" + 00 -1 

i=l -^0 ^ 



Notice that i?^, does not reduce to the standard 6' interaction on the line, see, e.g., 0, when it is restricted to 
a two-edge graph. Here we are following the notation in [22], |23j . The present 5' vertex is called sometimes 
5'g graph, where s is for symmetric. A discussion of the correct extension of the usual 6' interaction is given 
in [15j and |i8j. For completeness, we give the operator domain (the action is the same as in the other cases). 
We use the denomination 5' to avoid confusion with the previously defined interaction. 

ViH^^,) := U e H\g) s.t. ^^;.(0) = 0, ^,(0)-Vfe(0) = ^(VKO)-V'I.(0)) , j,k = 1,2, ...,n\ . 

Throughout the paper we restrict to the case of repulsive delta and delta-prime (5') interaction, i.e., 
a, /3 > 0. It is easily proved, for example by inspection of the operator resolvents in the following subsection, 
that such a condition prevents the corresponding Hamiltonian operator from possessing bound states. 

We also point out that the Hamiltonian Hf is well defined for a = and that indeed Hf\^_Q = Hp. We 

finally point out that, fixed a, ^ > 0, we shall consider the Hamiltonian operators Hf^ and H^/^ , namely, 
in the following we rescale a = av and (3 = /3/v. 



2.2 Resolvents, propagators and scattering 

For any complex number k with Im k > we denote by Rpik), R^{k), R^, the resolvents {Hp — k"^)^^ , {Hf — 
fe^)~^, {H^i — k'^)~^, respectively. 
We define the function Ut by 

e 4t 

Ut{x) := — = , t^O. 
In the following we shall use the same symbol Ut to denote the operator in L^(]R) defined by 

/oo 
Ut{x - y)i,{v)dy , 

Moreover, we define two integral operators acting on L?'{W^) 

r+oo 

ut : l2(M+) ^ l2(M+), [C/±V] (x) = / Ut{x ± y)^{y) dy , t ^ . 
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We stress that, according to our definitions, the operators and Ut have the same integral kernel, but act 
on different Hilbert spaces. 

For the cases we consider, resolvents and propagators can be easily computed (see, e.g. ^ pp. 201-226] 
for resolvent formulas with generic boundary conditions in the vertex). The results are summarized in the 
following theorem. 

Theorem 2.1. For any complex number k with Im A; > 0, the integral kernel of the resolvent operators 
Rpik), Rf{k), R^,{k) are given by 



RFik;x,y) 
R'iik;x,y) 
Ri{k;x,y) 



I 

2k^ 



2k 



■ ■ ik{x+y) 



2k 



2k a — 3ik 



_pi-k\x-y\j 

2k 2k3-iBk 





2 2 \ 


2 - 


-12 , 




2 -I J 


a — ik 


2ik 2ik 


2ik 


a — ik 2ik 


2ik 


2ik a — ik 


-1 + i 


Bk 2 


2 


-l + il3k 


2 


2 



2 
2 

■l + i^k 



(2.5) 
(2.6) 
(2.7) 



Furthermore, the unitary group of the related time evolution operators reads 

-12 2 



-iHpt 



U^l + Ut\\ 2 -1 2 
* 2 2-1 

2 
3 
2 



Ul,{x,y) = [Ut{x-y)-Ut{x + y)]l 
U^,^,{x,y) = [Ut{x-y) + Ut{x + y)]I 
where x,y G M+, a, /3 > 0. 



Ut{x + y) 



- / due-l^'Utix + y + u) 
3 Jo 



No 



+00 



due '3"f/t(j; + y + ii)JI, 



(2.8) 

J, (2.9) 
(2.10) 



Proof. We start from the proof of (2.6). Let = ipi for i = 1, 2, 3 and define R^ by 

f + 00 



iR^\ (x) 



2k 



M-y\i,,{y)dy + ^e'^-c, i = l,2,3 



where q = Ci{'^) are constant to be specified. It is obvious that (-R^)j satisfies 

■^-^')(m = i = 1,2,3 



(2.11) 



(2.12) 



then it is sufficient to fix a such that R^ belongs to T>{H^) in order to compute R^. The boundary 
conditions in (2.3) translate to the following linear system for Cj. 



e^'^^My) dy + ci 



+00 



z'^yMy)dy + 
^"'yMy)dy + 



C2 



C3 



(2.13) 



e'^^iMv) + My) + My)) dy - (ci + C2 + eg) 



JO 



+00 



e'^^My) dy + C2 
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the solution is easily computed and it is given by 



ci 



C2 



C3 



+00 



+00 



Jky 



k + ia 



3k + ia 
2k 



i'liy) + 



2k 



3k + ia 



2k 



3k + ia 



V'i(y) + 



3k + ia 
k + ia 



V'2(y) + 



2k 



3k + ia 

2k 
3k + ia 



V'2(y) + 



3k + ia 
2k 



V'2(y) 



3k + ia 
k + ia 
3k + ia 



i^3{y)j dy 
-03 (y) I dy 



(2.14) 



which gives (2.6). In order to obtain (2.5) it is sufficient to put a = in (2.6). Formula (2.7) can be proved 
by the same method. 



Now we prove (2.9). We start from the standard formula 

f+00 



UU^,y) 



1 
ni 



Rs{k; X, y) kdk, 



use (2.6) and recall the following identity 

1 gik{x+y) 

27r /_„ a — ik 



e-*'^'* dk 



+00 



(x+y+u) 

e 4t 
V'iirit 



(2.15) 



then we immediately arrive at (2.9). Formula (2.8) can be obtained by putting a = into (2.9). Formula 
(2.10) can be proved in the same way. □ 

Corollary 2.2. From the expression of the resolvent one immediately has the reflection and transmission 
coefficients: 



r{k) 
ruf (k) 
r^0 (k) 



1 

"3' 
k + ia 

3k + ia^ 

f3k + i 

I3k + 3i 



2 

~ 3' 
tH^{k) 



2k 



3k + ia 
2i 

' Bk + 3i' 



(2.16) 
(2.17) 
(2.18) 



We refer to |2T| for a comprehensive analysis of the scattering on star graphs. Indeed by using the results 
in [21] the reflection and transmission coefficients can be obtained directly by the boundary conditions in 
the vertex. 

Throughout the paper we shall need some auxiliary dynamics to be compared with the dynamics de- 
scribed by ( 1.1 ), so, for later convenience, we introduce the two-edge Hamiltonians Hj and the corresponding 
two-edge propagators e~*^^*, j = 1, 2, 3. 
Let Hj be defined by: 

V{Hj) := G H\g) s.t. ^,(0) = Vi+i(0) , ib'JO) + V'+i(0) = , ^^^(0) = , / j, j + 1} 




where in equation (2.2) it is understood that j = {1,2,3} modulo 3. 

The Hamiltonian Hj couples the edges j and j + 1 with a Kirchhoff boundary condition and sets a Dirichlet 
boundary condition for the remaining edge, so that there is free propagation between the edges j and j + 1 
and no propagation between them and the edge j + 2. 
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62 



■<- 

ei 




^i(O) = V2(0) 
^'i(0) + V2(0) =0 



V'3(o) = 



V>2(0) = V3(0) 

V^(o) + V3(o) = o 



^i(O) = 



i>3{0) = ^l(O) 

^^(o) + v'i(o) = o 



MO) = 



Figure 1: The open graphs depict the Hamiltonians Hi, H2 and ^^3. Under each graph we report the 
boundary conditions on vectors in the domain of the corresponding Hamiltonian. 



With a straightforward computation we have 



-iH,t 



where Tj are the matrices 



^0 10 
1 
.0 -1, 



To 



-1 0^ 
1 
1 0> 



^0 1^ 
0-10 
.1 0, 



(2.19) 



2.3 Strichartz's estimates 

A key tool in our method is the extension of the standard Strichartz's estimates (see e.g. |13j) to the 
dynamics on Q described by the propagators e~*^^*, e~^^s^, and e~^^s'^ . 

In this subsection we use the symbol H to denote any of the three Hamiltonians of interest described in 
section 12.11 



As a preliminary step, we remark that from equations (2.8), (2.9), (2.10), the standard dispersive estimate 
immediately follows: 

||e-^^**||^^^^||vl/||^„ t/0. (2.20) 

Proposition 2.3 (Strichartz Estimates for e~*^*). 
Let G L'^iG), r G L\L^, with 1 ^ g. A; ^ 2, | + ^ = |, and define 

^(t) = e-'^*^o Ht)= f dse-'"^'-'^T{s). 

Jo 

The following estimates hold true: 



for any pair of indices {r,p) satisfying 



I^IIlSL'- ^ cll^-ol 



\^\\ TP Tr ^ C r tQ Tk 



2 1 1 

2^p, r^oo, — I — = -. 

p r 2 



(2.21) 
(2.22) 

(2.23) 



The constants c in (2.21) and (2.22) are independent of T . 
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Proof. The proof is standard due to the dispersive estimate (2.20), see for instance [T3] and [19]. 



□ 



Remark 2.4. If H = Hf (H^,), the constants appearing in (2.20), (2.21) and (2.22) are independent of a 



(/5). Indeed, by the change of variable u 



ua [u 



u/P) the integral term in (2.9) ((2.10)) can be easily 



estimated independently of a obtaining a dispersive estimate (2.20) independent of the parameters and 
therefore, by the standard Strichartz machinery, uniform inequalities (2.21) and (2.22). 



3 Well-posedness and conservation laws 

Here we treat the problem of the well-posedness (in the sense of, e.g., [I3]), i.e., the existence and uniqueness 



of the solution to equation (1.2) in the energy domain of the system. Such a domain turns out to coincide 



with the form domain of the linear part of equation (1.1). Throughout this section, such a linear part 



is denoted by H, and, according to the particular case under consideration, it can be understood as the 
Hamiltonian operator Hp, Hf, or H^,. Correspondingly, we denote the associated energy domain simply by 
T>{£). All of the following formulas can be specialized to the particular cases 'D{£f), P(<£'^), or 'D{£^,). 

Let us stress that throughout the paper we do not approximate the dynamics in H^, but rather in L^. 
Furthermore, local well-posedness in is ensured by Strichartz estimates (proposition 2.3), as is easily seen 



following the line exposed in [13], chapters 2 and 3. Nonetheless, we prefer to deal with functions in the 
energy domain, since they are physically more meaningful. 

We follow the traditional line of proving, first of all, local well-posedness, and then extending it to all 
times by means of a priori estimates provided by the conservation laws. 

For a more extended treatment of the analogous problem for a two-edge vertex (namely, the real line with 
a point interaction at the origin), see [3]. 



First, we endow the energy domain "D^S) with the H^-norm defined in ( 1.5 ). Second, we denote by ^{S)* 
the dual of T>{£), i.e., the set of the continuous linear functionals on T){£-). We denote the dual product of 
r G T>{£)* and ^ £ T^i^) by (Fj^*). In such a bracket we sometimes exchange the place of the factor in 
T>{£y with the place of the factor in T>{£): indeed, the duality product follows the same algebraic rules of 
the standard scalar product. 

As usual, one can extend the action of H to the space 'D{£), with values in !){£)*, by 

{H^l,^2) ■■= (i/5^1,i/5^2), 



where (•, •) denotes the standard scalar product in L?'{Q). 
Furthermore, for any ^ £ the identity 

— e"^^*^ = -iHe-'"^^ (3.1) 
at 

holds in T>{£)* too. To prove it, one can first test the functional on an element H in the operator 
domain T>{H), obtaining 

/J \ / JHit+h)-^ _ iHt'^\ 



Then, the result can be extended to H G T^{£) by a density argument. 



Besides, by (3.1), the differential version (1.1) of the Schrodinger equation holds in 'D{£)*. 

In order to prove a well-posedness result we need to generalize standard one-dimensional Gagliardo- 
Nirenberg estimates to graphs, i.e. 

w^Wlp cw^'wlf^wnif^ +oo^p^i, (3.2) 
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where the C > is a positive constant which depends on the index p only. The proof of (3.2) follows 



immediately from the analogous estimates for functions of the real line, considering that any function in 
i/^(M+) can be extended to an even function in i/^(M), and applying this reasoning to each component of 

Proposition 3.1 (Local well-posedness in !){£)). For any G T^{£), there exists T > such that the 

and the following 



equation ([OJ) has a unique solution ^ G C°([0, T), P(f )) n C^([0, T), 
Moreover, eq. (1.2) has a maximal solution vj/™**^ defined on an interval of the form [0,T* 
"blow-up alternative" holds: either T* = oo or 



lim 



I iTfiiiax II 
l^t \\V{£) 



where we denoted by ^f^^^ the function ^ 



+00, 

"^^^ evaluated at time t. 



Proof. We define the space X := L°°([0, T), endowed with the norm 
Given \E'o G we define the map G : X — > X as 



SUPt6[0,T) 11^ 



t\\V{£)- 



e-'^-^o + i 



Notice that the nonlinearity preserves the space T>{£). Indeed, since any component ipi of ^ belongs to 
ir^(M+), then belongs to L2(M+) too, and so the energy space for the delta-prime case is preserved. 

Furthermore, the product preserves the continuity at zero required by the Kirchhoff and the delta case. 



By estimates (3.2) one obtains 



3 



so 



\G^\\x ^ \\^o\\v{e) + C / \\^ 



\t>{S) ^ 



l*ob(^) + cr||$||;^. 



Analogously, given [ 



IIG$ 



GE\\x ^ CT{\ml + \\E\\l) 11$ -H| 



X 



(3.3) 



(3.4) 



We point out that the constant G appearing in (3.3) and \^A\ is independent of ^'q, $, and H. Now let 



us restrict the map G to elements $ such that \\^\\x ^ 2||^'o||x)(£). From (3.3) and (3.4), if T is chosen to 



be strictly less than (8C||^'o|||,(£-0 ^, then G is a contraction of the ball in X of radius 2||^'o||x'(£')) and so, 



by the contraction lemma, there exists a unique solution to (1.2) in the time interval [0,T). By a standard 
one-step boostrap argument one immediately has that the solution actually belongs to C''([0, T), 2?(<?)), and 
due to the validity of ( |1.1[ ) in the space T>{£)* we immediately have that the solution ^ actually belongs to 
G^{[Q,T),V{£)))r\G\%T),V{£Y). 

The proof of the existence of a maximal solution is standard, while the blow-up alternative is a conse- 
quence of the fact that, whenever the P(£^)-norm of the solution is finite, it is possible to extend it for a 
further time by the same contraction argument. □ 

The next step consists in the proof of the conservation laws. 

Proposition 3.2. For any solution ^ G G'^{[Q,T),V{£)) f] G'^{[Q,T),V{£Y) to the problem the 
following conservation laws hold at any time t: 

W^tW = ll^oll, ^(*t) = ^(^-0), 
where the symbol £ denotes the energy functional 



1 



2 

Here the functional £iin coincides with£p, £f or£^,, according to the case one considers. 
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Proof. The conservation of the L^-norm can be immediately obtained by the vahdity of equation (1.1) in 
the space T>{£): 

d 



d II T ii2 



2Re(^,-^, 



by the self-adjointness of H. In order to prove the conservation of the energy, first we notice that {^t, H^t) 
is differentiable as a function of time. Indeed, 



t+h) 



and then, passing to the limit /i — )• 0, 

d 



2Im(|4',|2*,,if*, 



where we used the self-adjointness of H and (|1.1|). Furthermore 

d 



dt 



d^ 
dt 



(3.5) 



(3.6) 



From (3.5) and (3.6) one then obtains 

d 



dt 



and the proposition is proved. 

Corollary 3.3. The solutions are globally defined in time. 



□ 



Proof. By estimate (3.2) with p = co and conservation of the L^-norm, there exists a constant M, that 

1, 



depends on ^'o only, such that 



/ 1|2 



Therefore a uniform (in t) bound on is obtained. As a consequence, one has that no blow-up in finite 

time can occur, and therefore, by the blow-up alternative, the solution is global in time. □ 



4 Main Result 

In this section we describe the asymptotic dynamics of a particular initial state, which resembles a soliton 
for the standard cubic NLS on the line. 

According to section [3j we use the symbol H to generically denote the linear part of the evolution, 
regardless of the fact that we are considering the Kirchhoff, delta, or delta-prime boundary conditions. 
When necessary, we will distinguish between the three of them. 

We use the notation 

(j){x) = -v/2 cosh^"*" x , xGM, 
and for any xq S M and u E M we define 

(t>xo,v{x, t) = 6*5^6"**"^ e*V(a; -xo-vt), xeR, teR. (4.1) 

The function (pxo,v represents a soliton for the cubic NLS on the line which at time t = is centered in 
x = xo and has velocity v. Therefore, (j)xo,v is the solution of the integral equation 

(l)xo,v{x,t)=[Ute'r^(.-xo)]{x)+i[ [?7t_s|0^o,4-, s)| V^o,^(-, s)] (x)ds , t£R. (4.2) 

Jo 
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Let X be the cut off function x ^ C°°(M"'"), x = 1 in (2, +00) and x = in (0, 1). For later use we define 
also x+ = X[o,+oo) and X- = X(-oo,o]) where X[a,fe] denotes the characteristic function of the interval [a, h]. 
Moreover let xq and v be two positive constants and < 5 < 1. 
We take as initial condition the following function 



X{x)e *2^(/)(2; — Xq) 



■^o{x) 








XQ V 



1-5 



(4.3) 



and we denote by H,t the solution of the equation 

JO 



(4.4) 



The choice of the vector ^'o is used to render the idea that the initial condition is a soliton centered away 
from the vertex and moving towards the vertex with velocity v. The cut off function x in formula (4.3) is 



aimed at setting ^'o in the domain of the Hamiltonian H, see section 2.1 
Let us set t2 ■= xq/v + and define the following functions: 



^H,t2ixi,X2,X3) 



(xi + Xo - Vt2) 
!>(X2 - Xo + Vt-i) 



(4.5) 



t//e"^'^*2g»|x2gii2^(2.2 +X0- Vt2) 

Vi//e-*T-*2e-^t''3git20(^3 -X0 + uta)/ 

f t/fe-^'r^^g-if ^igit20(a;i - xq + ^^^2)' 




\iHi 



*T*2 g«5^3git2f 



(X3 + Xo - Vt2) 



(4.6) 



(4.7) 



They represent solitons on the line multiplied by the scattering coefficients of the linear dynamics rn and 
tn, that, in the particular regime we consider, are defined as follows: 



rHp = -1/3, tHp = 2/3 

_ 1 + 2za ~ _ 2 

3 + 2ia ^6 3 + 2id (4.. 



t - 



/3 + 6i /3 + 6i 



Remark 4.1. Notice that the coefficients f// and in can be obtained by the scattering coefficients (2.17), 
(2.18), identifying k with f/2 and replacing a by av and /3 by /3/w. This is due to the fact that we implicitly 
considered a particle with mass equal to 1/2, therefore the momentum k is linked to the speed w by /c = w/2. 

For any t > t2 we define the vectors ^■^jj ^ as the evolution of ^>'^ with the nonlinear flow generated by 
Hj, i.e., they are solutions of the equation 

= e-'^^^'-'^'>^^jjt2+^ f dse-'"^'^'-'^\^^Hs\^^Hs J = 1,2,3. (4.9) 
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Remark 4.2. The vectors ^''tt . can be represented by 



Hj.{xi,X2,X3) 



i^t2^it2jref ^^^^ 



e 4 



ef ( 



-X2 





2 



C-i.(-^3)/ 



where, for any t ^ 0, the functions (j)^'^'^ and i;^^'' are the solutions to the following NLS on the line 







(4.10) 



(4.11) 



ktr. 



tH 



Ut{x-y)e'^y(t){v-v^-^)dy + i I ds 



-oo 
oo 



Ut{x-y)e'^y(t>{v-v^-^)dy + i I ds 



Ut-s{x-y)\<p7f{y)\''€'Ay)dy 
Ut-s{x-y)\<pt{y)?<t^'I{y)dy. 



(4.12) 
(4.13) 



Our main result is summarized in the following theorem: 



Theorem 4.3. Fixed q;,/3 > 0, let H he any of the self-adjoint operators Hp, Hf" , H^!'" acting on L?{Q), 
where Q is the three-edge star graph, and defined by (2.1), (2.2), (2.3), (2.4). Call the unique, global 



solution to the Cauchy problem (4.4) with initial data (4.3). 

Then, there exist r* > and T,,, > such that for < T < T^, one has 



-t 



^l\\Ll{g) 

,1-5 



^ Cv ^ 



(4.14) 



for any time t in the interval xq/v + v < t < xq/v -\- v + Tlnu. 

In (4.14), C is a positive constant independent oft andv, the functions are defined by formulas (1.4), 
t and f, given in (4.8), are the scattering coefficients associated to H. 

The proof of the theorem will be broken into three steps or equivalently we break the time evolution of 
^H,t into three phases. 



Remark 4.4. A further consequence of Theorem 4.3 as in the case of |17| . is the fact that fast solitons have 
reflection and transmission coefficients which, up to negligible corrections, coincide with the corresponding 
coefficients of the linear graph. For example, a definition of the transmission coefficients along the branches 
j = 2,3 could be given considering the ratio between the amount of mass on the j— edge and the total mass, 
in the limit t — )• oo: 

lim 



t-s>oo II ^fJ 



\t'{v)\ , 



where denotes the restriction of the solution to the j— edge. In our case we do not have at our disposal 
the rigorous asymptotics for t — t- oo; nevertheless, we can obtain a weaker result. We have the results of 
Theorem 1.1, which give, in the time interval t2 < t < t2 -\- Tlnv , the estimate 



|ti| ) 



for a certain o" > and where P is the scattering coefficient of the linear Hamiltonian which describes the 
vertex. So, in the limit of fast solitons, i.e. u — )• oo, one can assert that the ratio which defines the nonlinear 
scattering coefficient converges to the corresponding linear scattering coefficient. And analogously for the 
case of reflection coefficient r(v), i.e. j = 1, one has (for t > t2, as before) 



hm ^ 



This is true for every coupling between the ones considered, i.e. Kirchhoff, 5 or 5' . 
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5 Proof of theorem 



4.3 



In the proof we drop the subscript H. When convenient, we specify the particular Hamiltonian operator we 
refer to. 

For any t G M we introduce the sohton 



„2 



Then, the function $t satisfies the equation 

$j = e-'^i^^-o + i / dse-'"^^^-'\^s?^3, t^O. (5.2) 
Jo 

5.1 Phase 1 

We call "phase 1" the dynamics in the time interval (0, ti) with ti = — —v^^. In this interval we approximate 



the solution by the soliton (5.1). The content of this subsection is the estimate of the error due to such an 



approximation, that is contained in proposition |5.3[ Before proving it, we need two lemmas. 
Lemma 5.1. Given ^ ta ^ tb ^ ti, for the functions 

Ki{t,x) := / Ut-ta{x + y)e-''^^y<i)[y-xo)dy + i ds Ut^^ix + y)e-'^ye-''^ e''(p^{y - xq + vs) dy 

Jo Jta Jo 

oo rt roo 2 



K2{t,x) := / Ut^tAx + y)e'^y4>{y + xo)dy + i ds Ut-s{x + y)e'^y e ^ e'' (^^y + xq - vs) dy 

Jo Jta Jo 

(5.3) 

the following estimate holds: 

\m\x,a,,m^)^ Ce--o+-*', i = l,2. 



where Xi„,,(M+) := L^^,^]L\R+) n LI^,^^L'{R+). 

Proof. Let us start with Ki. Adding and subtracting a contribution to negative values of y one can write 

/oo rt f'OO 2 

Ut-tSx + y)e-''^y(l){y-xo)dy + i ds Ut^,{x + y)e-'^ye-'''^ e''(p^{y - xq + vs) dy 
-oo Jta J — OO 

^ rt ^ . „2 . 

Ut-ta{x + y)e~'^y(p{y - xo) dy - i ds Ut-s{x + y)e-''^ye~'-^^ e'-^cti^iy - xq + vs) dy 

-oo Jta J —OO 

2 

e«5^e"**^e**(/)(x + xq - vt) 

rt rO 2 



Ut-ta{x + y)e ' (l){y - xo) dy - i / ds / Ut-s{x + y)e *2^e ^ e''(t)\y - xq + vs) dy 

(5.4) 



where we used the integral equation (4.2). By a straightforward computation, the X(^^ti,(IK''')-iiorm of the 
first term can be bounded by Ce"^^"^""**-*. To evaluate the size of the second term, let us write it as follows: 

oo 

Ut-ta{x-y)e'^^y<P{y + xo)dy =[Ut_t/'^-(t){- + XQ)]{x)= [Ut-taX+e'^^' <P{- + xq)]{x) , x > 0. (5.5) 
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Using the one-dimensional homogeneous Strichartz's estimates for Ut-tai namely, the analogous of (2.21) 
for functions of the half line, we can estimate the Xf^^tj^(M"'")-norm of this term as 

\\Ut-t^X+e'^^-H- + xo)\\x,^,^(R)^C\\x+e'r^(. + xon ^ Ce"^", 

where we used the notation Xt^,tjM) := L'^^^^^^L^{R) n Lf^^^^^^L^{R). 

The norm of the last term on the r.h.s. of equation (5.4) can be estimated in a similar way by 



ta 



[U.-s\4'-xo,v{s)\'^(l)-xo,vis)]ds 



V 



(5.6) 



Therefore, from (5.4), (5.5), and (5.6) we get 

\\K 



To estimate K2, the first term in its definition (5.3) can be treated as in (5.5), while the second is 
estimated following the line of (5.6). □ 

Lemma 5.2. Given ^ ta ^ tb ^ ti, let a and b two strictly positive numbers, with 



b ^ 



1 



(5.7) 



8a2 + 4a ' 

Moreover, let y be a real, continuous function such that ^ y{ta) ^ cl, o,nd 

^ y{t) ^ a + by'^{t) + by^{t), for any t G [taM- 
Then, 

max y{t) ^ 2a. 

t(i[ta,tb] 

Proof. Consider the function fb{x) := bx^ + bx"^ — x + a. Denoted b := ^^2^4^^ , one has fi{2a) = 0. 
If 6 ^ 6, then fb{x) ^ fii^) fo^ ^'^Y x > Q. Besides, notice that /fe(0) = a > 0, then there must be a point 
X G (0,2a] s.t. fb{x) = 0. Finally, since the function y is continuous, in order to satisfy the constraint (5.7) 
one must have 

y{t) ^ X ^ 2a, for any t G [ia^^ft]- 

□ 



Proposition 5.3. Let be the solution of the equation (4.4), and (^t be the solution of equation (5.2). 
There exists C > 0, independent oft and v, such that 



\'^t-<^A\ ^ Ce" 



(5.8) 



for any t £ [0, ti]. 



Proof. Let us define := — ^t, and fix ta G [0,ti]. Then, from equations (4.4) and (5.2), we have 



-iH{t-ta): 



ta + (e"'^(*-*") - e-^^i^*-*"))-^*, + i Ae-^-^^*-') - e-'^i(*-"))|$,p$, ds 

Jta 



+ i f e-^^(*-^) [iH.pH, + + |$,|2h, + 2Re(H,$,)H, + 2Re(H,$,)$, 

Jta 

e-^H(^-'^)Et^+F{ta,t) 



+ i 



-iH{t-s) 



+ + + 2Re(H,$,)H, + 2Re(H,$,)$J , 
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where we defined 

Let us fix tb £ [ta,h], and denote Xt^^tt = -^[^ tt]-^'^ ^ ^fta tb]-^^' '^'^^^ 

,-iH(-s) r|s^|2s^ + |2^|2^^ ^ l^^p^^ ^ 2Re(H;$,)H, + 2Re(H;$,)$, 



ta 



(5.9) 



Using (2.21) the first term on the r.h.s can be estimated as 



the cubic term. Since both pairs of indices (oo, 2) and (6, 6) fulfil (2.23), in (2.22) we can choose q = k = 6/5 
and obtain 



We estimate the integral term on the r.h.s. of (5.9) also using Strichartz's estimates. Let us analyse in detail 

= 6/5 

(5.10) 



Jta 



^ C|||H.| H. 11^6/5 



Moreover, by standard Holder estimates, 



|2 I 

Il6| 



l^l|2 



Then, by (5.10) and (5.11) 



ta 



(5.11) 



(5.12) 



Notice that the constant C can be chosen independently of ta,th, and of the boundary condition at the 
vertex. 



The other terms in the integral on the r.h.s. of (5.9) can be estimated analogously. One ends up with 

mUa,, ^ c iiHtji + \\F{ta, oiix.^,, + c{h - tS^mwi^^^^^ + ^(^ - tjmwi^^^^^ + cit^ - tjmwx,^^,^, 

where the arising norms of <I> were absorbed in the constant C. 

If th and ta are sufficiently close, then C{th — ta)^ < 1/2. Furthermore, since the quantity th — ta is upper 

2 1 — 

bounded, one can estimate (4 — ta)^ by C{tb — ta)^ ■ So, for some C > 



^ C \\EtJ+C\\F{ta,-)\\x,^,+C{t,-tar^ 



+ 



(5.13) 



Applying lemma 5.2 to the function y{t) = \\E\\xt^ t, which is continuous and monotone, one has that, if 

-2 



tb-ta ^ [SC ( \\EtJ + \\F{ta,-)\\x,^J' + 4C' imtJ + \\F{tar)\\x,^J 

then ||S||xf^,tj^ ^ 2C \\Et^\\ + 2C \\F{ta.,-)\\xta,%- From the immediate estimates 

\\EtJ ^ 4, ||F(t„.)|k,,,, ^ limOIko,,, 

if one denotes 

r := fsC^ (4+||F(t„.)|Uo,j' + 4C2 (4+||F(t„,.)||xo.J 
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then for any te [0,ii) ^ 2C \\'Et\\+2C 

We divide the interval [0, ti] in + 1 subintervals as follows 



[0,ti] = (uf^-o'b'^'(j + l)^])u[iVr,ti], 



where 



N :-- 



Xq-V 



1-5 



VT 



integer part. 



Making use of lemma 5.2, and noting that ^ (j^^)^, one proves by induction that 

^ (2C^y+'Po|| + E(2^)^'^'"'ll^(^^'-)lk^^^^^ j = 0,...,iV-l, 



k=0 
N-1 



(5.14) 



^ (2C)^+iSo|| + (2C)^+i-'=||F(fcr, OIU,.,,^,,. + 2C \\F{Nr, .)\\x 



NT,ti ' 



k=0 



where the last inequality comes from the fact that ti — Nt ^ r, so lemma 5.2 applies to this last step too. 
The norm of S as a function of the whole time interval [0, ti] can be estimated by 



N-1 

I'^II^O.ti ^ ^ ^ + ||"| 

j=0 

3=0 j=0 k=Q 



(5.15) 



In order to prove the theorem using (5.15), we need more precise estimates for ||Ho|| and [[^(jr, •)||xj^ {i+i)T- 
First, 

-2X0 



l^of ^ / (p'^ix- xo)dx+ / (/)2(x + xo)= 2(1 -tanh(xo - 2)) ^ Ce~ 
Jo Jo 



(5.16) 



To estimate ||-F(ia, Oll^t^.t specialize F to the three cases under analysis. From the explicit propagators 
d^Igl), (|2T0l), we get' 



Fl/^{ta,t) =Utt^ 1-1 1 Q I U^ 



-1 -1 2 



FF{ta,t) =^Ul^^ I -1 -1 2 I ^t^ + 

2 2 2 

Ft{ta,t) = Fp{ta,t) - -av / dne-t™ {U.ttJ^t.) (• + 



-1 -1 2 

-1-12) \<^s\^^sds 
2 2 2 



2 

i-av 
9 



+ 00 



ds due-t™ ([/+^JI|^.,|2$^) (. + n) 

ta Jo 



1 -1 



ta 



2v 
'J Jo 
It is immediately seen that 



2 

+00 _ 3 

due r''{U+tJ^ta)i- + 



1 -1 

-1 1 I \^s\^^sds 
2 



. 2v 
u) — I— 



I ds due r^{U+J\<^s\'^s)i- + 

P Jta Jo 



u 



FF{ta,t,x) 



-Kiix,t) - K2ix,t) 
-Ki{x,t) - K2{x,t) 
2Ki{x,t) + 2K2{x,t), 
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where Ki and K2 were defined in (5.3). Lemma 5.1 yields 
Furthermore, since 



(5.17) 



2 ^+00 _ / Ki{x + u,t) + K2{x + u,t) 
F^''{ta,t,x) = FF{ta,t,x) - -dv due't™ Ki{x + u,t) + K2{x + u,t) |, 

\ Ki{x + u,t) + K2ix + u,t) 

after the change of variable u ^ uv we conclude 



(5.18) 



Finally, 



Fg/ (ta,t,x) 



Kiix,t) - K2ix,t) 
-Ki{x,t) + K2{x,t) 




2y r+oo _3^^ [ Ki{x + U,t)+K2{x + U,t) 

due /3™ Ki{x + u,t) + K2{x + u,t) 
^ -^0 V Kiix + u,t) + K2{x + u,t) 



yields, after the change of variable u — )• uv, 



(5.19) 



Now we go back to estimate (5.15 ). Due to (5.16), (5.17), (5.18 ), and (5.19 ), and estimating any geometric 



sum by the double of its largest term, which is justified if the rate of the sum is not less than two, we get 



i=o j=o fe=o J 



N-l 



k=o ^2Cj 



■xo+vti 



(5.20) 



N~l 



^ 2C(2C')^+^e-^« + 4CCe-^o+^^ ^ e^'"^ + 8CC^e"^o+^^^ + 2CCe-^"+''*i 

j=0 

^ 2C7(2C')^+^e-^« + 8CCe-^o+^''^ + 8CC2e-^°+^''^ + 2C7Ce~^°+^*^ 



Concerning the first term on the r.h.s. of (5.20), we have 



(2C7)^e-^o ^ 



(5.21) 



From (5.20) and (5.21) we get ||H||xoi^ ^ ^ , so (5.8) follows and the proof is concluded 



□ 



5.2 Phase 2 

We call "phase 2" the evolution of the system in the time interval (ti,t2) with t2 = ^ + v~^. Let us define 
the vector 
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Figure 2: We represent the state for t = ti and t = t2. The continuous Unes stand for the real graph, 
the dashed hues represent the extension of each edge to negative values of x. At time t = ti only the vector 
^s,in -g gignificantly supported on the real graph. At time t = t2 only the state $■5'.°"* j^as relevant support 
on the real graph, while the body of the soliton in has moved to the negative axis associated with the 
first edge. 



with 



S.in 



((kxo-v{t)\ 



\ / 



S.out 



(r(t>-xo,v{t)\ 
i4>-xo,v{t) 



where the function cl)xo,v was defined in equation (4.1 ) and the reflection and transmission coefficients, f and 



t, must be chosen accordingly to the Hamiltonian H taken in the equation (4.4). The explicit expressions 



of f and t in all the cases H = Hp, Hf", H^/'" can be read in formula (|4.8l). 



Proposition 5.4. Let t G (^1,^2) then there exists vq > such that for all v > vq 

ll^'t - $f II ^ Civ~^ , 



moreover 



s.out I 



't2 - ^t2 II ^ ^^2^ 2 , 

where Ci and C2 are positive constants which do not depend on t and v. 



Proof. From the definition of see equation (4.4), we have 



e-*^(*-*i)^t, +i f dse-'^^^-'^l^sl'^'^s- 



We start with the trivial estimate 



(5.22) 
(5.23) 



(5.24) 



and estimate the r.h.s. term by term. The estimates involved in the analysis of the first term are similar to 
the ones used in the previous proposition, thus we omit the details. Similarly to what was done above we 
set Xtj^^t2 = -^^[^j 12]^'^ ^ ^fti t2]'^^' ^^^^ Strichartz estimates (see equation (5.12) and proposition 



2.3) 



dse-'^^-'M^J^^. 



^ C{t2-ti)^/^\M 



X, 



tl.t2 



Xti,t 



(5.25) 
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which imply 
By lemma 



5.2 



one has that if (^2-^1) ^ [SC^H^rtJp + 4C2||^'jJ|] , then \\^\\x,^^.^ ^ 2C||^'ti ||; using this 

(5.26) 



estimate in the inequality (5.25) we get 



II II 11^*1,(2 II II 



where we used t2 — h = 2v 



We proceed now with the estimate of the second term on the r.h.s. of inequality (5.24). Let us set 

4'-xo,v{tl) 

\ J 



^1 



We notice that $ 



adding and subtracting $ 



where the vector $t was defined in equation (5.1) and rewrite by 



The following trivial inequality holds true 



\^-iH{t-t^)^ _ ^-iH(t-ti)^S,in\\ ^ 



(5.27) 



5.3 



and the fact that ||^>^f'|| ^ 2^" 



where in the latter estimate we used proposition 

Let us consider now the last term on the r.h.s. of inequality (5.24). We are going to prove that for all 
t £ (^1,^2) and for v big enough 



Let us introduce the functions 



and 



/•oo 

4>ti^)-= / Ut^td^ - y)4'xo-v{y,ti)dy 
Jo 



(5.29) 



(5.30) 



First we prove a preliminary formula for the vector e ti)^^'*" ^ggg equations (5.35) and ( 5.36| ) below). 
For any constant a > 0, not dependent on v let us consider the term 



va 



-^^''[Ut.t,^,,,^v{ti)]{u + x)du 

/•CO 

--va I du dy e-™''[/t_ti {u + x + y)e*'^(*i)e-^t?'(/,(y -xo + vh) 
Jo Jo 

2 

where we set ip{t) := —t^ + 1. By integrating by parts we obtain the equaUty 

"'[Uttt,(^.o,-vih)]iu + x)du 



va 



4 

oo 



--2ia / du dye-"'"'[/t_t,(u + x + y)e'^(*i) 
Jo Jo 

--Ai{x, t) + A2ix, t) + A3{x, t) 



dy 



-e 2 



<j){y -xo-\- vti) 



(5.31) 
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with 



/oo 



A2{x,t) := -2ia 

A-i{x,t) := -2ia 
We notice that 



d 



— Ut-tA'^ + x + y) 



du I dye^ 

Jo L"y 

/OO 
dy e-'''"'Ut-t, iu + x + y)e^^(*i)e-*5S/ 



due-"™ + 

' L2(R+) 

Then the foUowing estimate for the term Ai holds true 



due-'''"'Ut-t,{u + -) 



dy 

Wx+e 



(t>{y -xq + vti) 



1 

2va 



Pi(i)llL2(R+) ^2a 
The term A3 is estimated by 



,,1-5 



\(^{-xo + vti)\ ^C—j^. 



(5.32) 



\\Mt)\\mR+) ^2a / du 



-uva I 



^2 



Il2{r) 



a I due 



[Ut-tiX+e '2'</''(- - xo + vti)] ( - (n + •)) |L2(]r+) 
C7 



(5.33) 



where we used the equahty [Uf^ f]{x) = [UtX+f]{—x). 
We compute finahy the term A2. By integration by parts 



A2{x,t) 



2ia / du 







dye' 



-Ut-tAn + x + y) 



ei'/'(*i)e-*!?^(/,(y - xo + uti) 



■2iact)t {x) - 2iva^ / e-'^''"'[U^_^^(t)^^^.^,{tl)]{u + x)du , 
Jo 



where the function (pf was defined in equation (5.30). Using the last equality in the equation (5.31) we get 



va 



' [u+_t^^^,,.^,{h)] (m + x)du = Y^r^^tix) 



2ia . Ai{x,t) + A'i{x,t) 



+ 



1 + 2ia 



(5.34) 



From the definition of f^f'^"' and using the last equality with a = a/3 in the formula for the integral kernel 

TTOtVJ- J . I . „ 



of e see equation (2.9), it follows that 



2 - + 



V 



Aa{t) 

\Mt)J 



(5.35) 



V / 

where the function (pj' was defined in equation (5.29) and we set Aa{x,t) := [{Ai{x,t) + A3{x,t))/{1 + 
2^0)] Similarly using equality (5.34) with a = 3//3 in the formula for the integral kernel of the 



propagator e , see equation (2.10), we get 













^-iHl/\t-H)^S,in ^ 


l3+6i ^* 




2 

~ 3 








J 







(5.36) 
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. We notice that the 



where we introduced the notation := + A^{x,t)) /{I + 2m)] j^^g^^ 

analogous formula for e'^^P^^-^^^^^f" can be obtained from equation (|5.35p by setting q = and = 0. 

To get the estimate (5.28) we show that, at the cost of an error of the order (t2 — ii), for t S {ti,t2), the 
functions 4>~[{x) and cf^tix) can be approximated by the solitons (t)xo-v{x,t) and 4>-XQ,v{x,t) respectively. 
We consider first the function (pf , by adding and subtracting a suitable term to the r.h.s. of equation (5.30) 
we get 



Ut-tiix + y)(l)cco~viy,ti) + i ds I Ut-s{x + y)\4> iy,s)dy 



Ut-t^{x + y)(l)xa,^v{y,ti)dy - i j ds j Ut-s{x + y)\(l)xQ-v{y, s)\'^4>xo-viy, s) dy 

•j Jtl J — oo 

=(j)-xo,v{x, t) + I{x, t) + II{x, t) , 

where WG usGci the fact tha-t (px^^—v 

{—x,t) = (j)-xo,v{x,t) and we set 



(5.37) 



I{x,t) :-- 



Ut-ti {x + y)(f)xo-viy, ti)dy 



and 



II{x,t):=-i ds Ut-six + y)\(j)xo-v{y,s)\'^(l)xo-viy,s) 

Jtl J —oo 



dy. 



For the term I we use the estimate 



I-^IIl2(r+) ^ 



The term // is estimated by 



xo + vti)\\i2^^-) ^2e " 



Similarly, for the function (j)^^ , we get 

(pJix) = ^xo,-v{x,t) + III{x,t) + IV{x,t) 

where we set 



(5.38) 



III{x,t) := - / Ut-tiix - y)4>xo-v{ti,y)dy 

J — oo 



and 

/^t /'OO 

IV{x,t):=-i ds Ut-six - y)\(l)xo-vis,y)\'^(l)xo-vis,y)dy . 

Jtl J —oo 

For t G (ti,t2) the estimates 



\IV\ 



L2(K+) 



^ Cv 



-5 



are similar to the ones given above for the terms / and //. Then from equations (5.37) and (5.38) we have 



xo-v{l^)\\L2{R+) 



for all t G {ti,t2)- By using the last estimates and the estimates (5.32) and (5.33) in equations (5.35) and 
(5.36) we get 



^^_iH(t-ti)^S,in _ ^ 



yl/2 y 



which in turn implies that for v big enough the estimate (5.28) holds true. 



Using estimates (5.26), (5.27) and (5.28) in the inequality (5.24) we get the estimate (5.22). 
The estimate (5.23) is a consequence of estimate (5.22) and of the fact that jj^f^ ~ ^t2 



S.oui I 



1$ 



S,in I 
t2 I 



Ce'" 



□ 
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Remark 5.5. Notice that estimate (5.22), although not strictly necessary for the proof of theorem 4.3 



enforces the picture that in the phase 2 a scattering event is occurring. The true wavefunction can be 
approximated by the superposition of an incoming and an outgoing wavefunction. At the end of this phase, 
only the outgoing wavefunction is not negligibile. 

5.3 Phase 3 

Let us put = t2 + Tin 11. We call "phase 3" the evolution of the system in the time interval {t2,t3). The 



approximation of during this time interval is the content of theorem 4.3 
We recall the following result (see [17J). 

Proposition 5.6. Let (f>l^ and (t)^^ he defined as it was done in equations (4.12) and (4.13) above. Then 
Vfc G N there exist two constants c{k) > and (T{k) > such that 

c{k)(lnvy(''^ 



^^t\\L^iR-) + \\4>t\\L° 



(5.39) 



for 7 = {ref, tr}, uniformly in t G [0,Tlnf]. 

Notice that also the norms ||</'7llip(iR-)) for 2 ^ p ^ oo, are estimated by the r.h.s. of (5.39). Now we 



can prove theorem 4.3 



Proof of Theorem 4-3. The strategy of the proof closely follows proposition 5.3 We will just sketch the 
common part of the proof while proving in details the different estimates. 



Let us define := — X^ i=i '^l where the vectors <I>j were given in equation (4.9) and fix ta G [i2ii3]- 



From equations (4.4) and (4.9) it follows that the vector satisfies the following integral equation 



7, -p-iH{t^ta)'^ 



ds ( e 



+ ^ (^e-^^(*-*«) - e-^^^ <1>^'^ + i I 

j = l L Jta 



'iH{t-s) 



Jl J2 J3 



+ i I dse-*-^(*-') 

'ta 



3 2 3 r 3 n 3 

j=l j=i L J 

e-'^(^~^^)Et^+G{ta,t) 

+ i rdse-^(*-)[H, + ^cI>i'H, + |H,|2^$^, + 2Re ^ ci>^,] ^ ci>^^ 

tn. _ ^ 1 1 _ ,■ 1 _ ,■ 1 



(5.40) 



i=i 



Let us fix tb £ [tajts] and let Xt^^t,, = L'^ M ^"^ ^ ^ft t^]^^' Using the Strichartz estimates as it was done in 



proposition 5.3 (see equations (5.9) - (5.13)), it is straightforward to prove that 



^ c 



'^u\\ + \\G{ta,-)\\x,^,,+{h-taf'^ 



+ 



|3 



where C depends only on the constants appearing in the Strichartz estimates. Using lemma 



5.2 



as it was 



done in the proof of proposition 5.3 it follows that there exists r > such that, for any t £ [^25^3) one has 

^ 2C \\Et\\+2C\\G{t,.)\\x,,^^. 
We divide the interval [^2,^3] in + 1 subintervals: [^2 + JT, ^2 + (j + 1)''")) with j = 0, ...,N — 1; and 
[t2 + Nt, t^), and where N is the integer part of (^3 — t2)/''"- 
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Then proceeding by induction as we did in the proof of proposition 5.3, see equations (5.14) and (5.15), we 
get the inequahty 



N N j 

j=0 j=0 k=0 



(5.41) 



Now we estimate the initial data ||Ht2|| and the source term ||G'(ta, Oll^ta with t2 ^ ta ^ tb ^ and 
tb-ta^ T. 



By proposition |5.4| (estimate (5.33)) and using the definitions (4.5) - ( |4.7[ ) one has 

/'i(j)xo-v{t2)\ 
f(t)xo-v{t2) 



-t2 



\i(j)xo-v{t2)J 



with 11^*2 - ^fj""*!! ^ Cv-^/^. Since 



\(pxo-vix,t2)\'^dx 



{x — xq + vt2)\'^dx = 2 



roc 

/ sech( 



x)'^dx 



4e- 



1 + e-^^^ 



^ 4e' 



we have 



(5.42) 



Let us now consider the source term G{ta,t). We use the estimate ||G'(fa5 •)||xt^,tj^ ^ l|G'(t2) •)l|xt2,t3 • To 
simphfy the notation we set G{t) = G{t2,t) and 

.■1 L ^ to 



G2{t):=i d.e-^(*-^) (^-^nn^hjs)^^^^^^^^- 

J t2 



Jl J2 J3 



By the definition of G{ta,t), see equation (5.40) it follows that 

G{t) = Gi{t) + G2{t) ; 

we estimate Gi{t) and G2{t) separately. 

We proceed first with the estimate of the term Gi. From equations (2.8), (2.9), (2.10) and (2.19) one 
can see that for any (column) vector F = {Fi, F2, F3) G 

(U^Fi\ (va e-™'* Pt ^i] + ■)dy\ 

2 



-iHt „-iHA 



e-'"^^F 



UtF2 
WF3J 



va 



/g°°e-"™[C/+F2] {u + -)du 
\va e-™'* [U^Fs] {u + ■)duj 



j = 1,2,3; 



where the constant a and the matrices M^- must be chosen accordingly to the Hamiltonian H: 



for H = H^", 



Ij = -I + -J - Tj ; 
o 
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for H = hI/\ 



and the formula for H = Hp can be obtained by setting d = in the formula for H = Hf 
Then, denoting by / = 1,2,3, the /-th component of the vector Gi one has 



with 



(Gi(x,t)),:= (M,)^, 
j,k=i 



t2 



2va 



E 

j,k=i 



[C/,t,,K)J(n + x) + . / ds[Ul,\[^i)^\\^i)^]{u + x) 



t2 



du . 



From the definition of the vectors see equations (4.10) - (4.11), we see that for each I = 1,2,3 the 
function is a linear combination of four functions, f^' and f^'~ , with 7 being equal to ref and 

tr, given by 



ft (x) :=e""4-e" 
and 



Ut-t,{x + y)(t>l{y)dy + i ds Ut-s{x + y)\4>]_t^{y)\'^^]_t,{y)dy 

Jt2 Jo 



(5.43) 



Ut-t-Ax + y)4>l{-y)dy + i ds Ut-s{x + y)\<fi_tA-y)\^(fl_tA-y)dy 

Jt2 Jo 



where the functions c/x^^^ and cfk'^'^ were defined in equations (4.12) and (4.13) respectively. 

Similarly one can see that for each / = 1, 2, 3 the function (Gi(x, t))^ is a linear combination of 



va 



va 



fl^^'~^ {u + x)du , va 



'Vua ftr,+ 



f^' {u + x)du 



fl^^'- (n + x)du , va I e-™"/f {u + x)du . 



First we study the function ff'~^ ■ We notice that, adding and subtracting a suitable term in equation (5.43) 
and using the definitions (4.12) and (4.13), ff'~^{x) can be written as 



= I{x, t) + II{x, t) + III{x, t) , 



with 



II{x,t) := -e-'^^^e'^^ 



III{x,t) := -ie-'^^^^e'^^ 



t 

ds 

t2 JO 



Ut^t2{x - y)<pQ{-y)dy 



Ut.s{x - y)\4>UA-y)\^UA-y)dy . 



Similarly to what was done above, we set Xt^^tii^"^) = L'^^ ^^jL^(M='=) n L^^^ (M='= ) . By proposition 
we have 

c'(/t)(lnv)'^'('=) 

/TnU-\ — lie/; , II •\z /jm^\ :\ - 



5.6 



\Xt. 



yk{l-5) 



(5.44) 
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where c'{k) and cr'{k) are constants, different from the one appearing in proposition 5.6 For our purposes 
we do not need to compute them. 

Using the one dimensional Strichartz estimates for Ut, we have 



\^I\\xt^,tJM.+) ^ C'IIx-'/'oIIl2(k) ^ Ce 



(5.45) 



Finally, the term /// can be estimated using the inhomogeneous Strichartz estimate and proposition 5.6 



,,fc(l-<5) 



(5.46) 



Collecting the estimates (5.44), (5.45) and (5.46), it follows that 



II Z-"^' IUt,,t3(R+) 



^ c 



The estimate of /7' is similar and we omit it. We have proved that for some d{k) and cr' {k) possibly bigger 
than c{k) and (j{k) we have: 

'd{k){\nvY'^^y 



|Gl||Xt2,i3 ^ C 



where Gi is the vector in with components (Gi(a;,t))p / = 1,2,3. 

The estimate of Gi{t) = ((Gi(t))-^, {Gi{t))^, (^i(*))3) ^ trivial consequence of the fact that 



va 



-vua f7,it 



/.T'±(u+ ■)du 



'7, it I 



X-, 



^ C 



z{k){lnvY(''y 



from which it follows that 



and 



\G-\ Wxt , ^ C 



c'{k){lnv 



d{k){\nvY'^^y 



\\Gi\\xt^,t., ^ C 

We analyse now the term G2. Due to the presence of (1 — 5jxj2^j2n) the components of the vector 
contains only terms (up to a phase) like 



(5.47) 



where 71, 72 and 73 can be ref or tr. This can easily be seen by using equations (4.10) - (4.11). By 
Strichartz methods, it is sufficient to estimate the Lj^^ j^jL^(M"*') norm of these terms. Then using Holder's 
inequality we have, for istance 

|iei,et.'^7^t,(-)llL^(R+) ^ lie.jL4(R+)|ietjL4(R+)lietjL°°(R-) ^ C ^ ^^^^^ . 

The second kind of terms can be estimated in the same way and we obtain 

'd{k){\nvY'^^y 



|G2||xt„.t, ^ C 



yk{l-5) 
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which, together with the estimate (5.47), gives 



Fix k such that k(l — 5) > 2 then for v sufficiently large (5.48) implies 



(5.48) 



(5.49) 



N N j 

j=0 j=0 k=0 

N 

^2C^;-V2(2C')^+i + 2v-' Y.^2Cy^^ 

j=0 

i^2Cv~^/\2C)^+' + 2v-\2C)''+^ ^ Cv-^'\2C)'' 
Since is the integer part of {t-^ — t2)/T = ^ Inv we have 



We can finally set r* = T/ln(2(7) and = (5r*/2 and obtain 



which concludes the proof of theorem (4.3). 



□ 



6 Conclusion and perspectives 

In the present paper we have given a first rigorous analysis of nonlinear Schrodinger propagation on graphs. 
We have given a preliminary proof of local and global well posedness of the dynamics, and of energy and 
mass conservation laws for some distinguished vertex couplings, i.e Kirchhoff, 5 and 5' couplings. Then 
we concentrated on the problem of collision of a fast solitary wave on the graph vertex (with couplings as 
before) . It turns out that the solitary wave splits in reflected and transmitted components the form of which 
are again of solitary type, but with modified amplitudes controlled by scattering coefficient given by the 
linear graph dynamics. This behaviour holds true over times of the order In v where v is the velocity of the 
impinging soliton. 

We add some other remarks on the result and further analysis and generalizations. 
To begin with, let us note that the real line with a point interaction at can be interpreted as a degenerate 
graph with two edges. The cited paper [T7] treats the special case of a 5 interaction on the line, and our 
description shows how it could be possible to extend their results to other point interactions; among the 
examples treated in the present paper there is a version of the 5' interaction, showing how to treat point 
interactions of a more singular character than the one given by a (5. 

Concerning more general issues, a sharper description of the post interaction phase can be achieved by 
an explicit characterization of the evolution of the modified solitary profiles, i.e. of the . This last part is 
somewhat delicate, and intersects with contemporary intense work on asymptotics for solitons in integrable 
and quasi integrable PDE, so we limit ourselves to the following remarks. In the case analysed in [T7], the 
asymptotic behaviour of nonlinear Schrodinger evolution of solitary waveforms with modified amplitudes 
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is given, and making use of inverse scattering theory it is shown (appendix B of the cited paper) that the 
evolution is close to a soliton up to times of order In v and an error the order of which is an inverse power of 
V. Borrowing from these results, it is possible to get in our case too, but we omit details, the asymptotics of 
the i.e. of the free NLS evolution of the modified solitary profiles outgoing from the phase two. It turns 
out that these outgoing wavefunctions can be approximated, on the same logarithmic timescale of Theorem 
1.1, as new solitons with the same waveform of the unperturbed dynamics, modified amplitudes and phases, 
plus a dispersive ("radiation") contribution. The meaning of this statement is that the L°° norm of the 
difference between the evolved modified solitary profiles $^ and such final outgoing solitons has the usual 
dispersive behaviour, \t — ^2!"^. Let us note that, following this strategy, at the end of the phase three, 
there would be two types of errors: errors due to the approximation procedure in phase one and two {O12); 
and errors arising from neglecting dispersion in the reconstruction of the outgoing solitons (Oioo) . 

An important question concerns the possibility of extending the timescale of validity of approximation 
by the solitary outgoing waves. As a quite generic remark, this possibility could be related to the asymptotic 
stability of the system, or of systems immediately related to it. 

More concretely, in a different type of model (scattering of two solitons on the line) in the already cited 
paper pQ, some considerations are given on obtaining longer timescales of quasiparticle approximation in 
dependence of the initial data and external potential, but it is unclear whether similar considerations can 
be applied to the present case. 

Another issue is the nonlinearity. The fundamental asymptotics proved in [TTj and used in the present paper 
relies on the integrable nature of cubic NLS, and it is not immediate to extend these results to more general 
nonlinearities. One can conjecture that for nonlinearities close to integrable ones which admit solitary waves, 
the outgoing waves are close to solitons over suitable timescales. Let us mention, however, the recent results 
of Perelman on the asymptotics of colliding solitons for nonlinearity close to integrable or critical on the 
line ([52] and [23j ). 

A final problem is the extension of results of the present work to more general graphs. We believe 
that results similar to the ones of the present papers are valid for more general boundary conditions at the 
vertex of a star graphs, with the same proof, under the condition of absence of eigenvalues for the linear 
Hamiltonian describing the graph. In presence of eigenvalues, some Strichartz estimates weaken, and a more 
refined analysis is needed (see |14] for the analogous problem on the line with an attractive 6 interaction). 
Of course, the extension of the present results to the case of star graphs with more than three edges has 
to be considered straightforward, while the extension to graphs having a less trivial topology is an open 
problem. 

Finally let us comment briefly the recent paper [26 . In this partly heuristic paper the authors study 
a star graph (but also more general type of graphs are considered) with a NLS in which on every edge 
there is a different strenght /3k in front of the the cubic term. The authors fix a boundary condition which 
guarantees that mass and energy of the solution are (formally) conserved. Moreover according to the authors 
it is possible to derive a condition on the strenghts Pk which allow for complete transmission of an incoming 
solitary wave across the vertex. In these same situations the authors show that an infinite chain of conserved 
quantities exists, defined analogously to the case of the NLS on the line. The result, if formal, is interesting, 
and concerning the relation with ours we note the following. In the case of a three edge graph and more 
generally for a odd edge number, the complete transmission is made possible exactly by the fine tuning 
of the coupling constants in front of the nonlinearities. For the case of a single medium with the same 
nonlinearity on every edge and Kirchhoff boundary conditions, one can prove (see [2] , where more generally 
the case of nonlinear bound states for 6 boundary conditions is treated) that exact travelling solitons exist 
only in the case of a graph with an even number of edges, while in the case of an odd number of edges a 
stationary state is formed which is given by half a free soliton on every edge. 
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